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Abstract. In this paper we study algorithmic aspects of tropical intersection 
theory. We analyse how divisors and intersection products on tropical cycles 
can actually be computed using polyhedral geometry. The main focus of this 
paper is the study of moduli spaces, where the underlying combinatorics of 
the varieties involved allow a much more efficient way of computing certain 
tropical cycles. The algorithms discussed here have been implemented in an 
extension for polymake, a software for polyhedral computations. 



1. Introduction 

Tropical intersection theory has proved to be a powerful tool in tropical geome- 
try. The basic ideas for an intersection theory in R™ based on divisors of rational 



functions were first laid out in [Ml and were then further developed in AR . Even 
earlier, [FS2| had proved the fan displacement rule in the context of toric geome- 
try. It describes how cohomology classes of a toric variety -^(A) can be multiplied 
using a generic translation of A. This was later translated to the concept of stable 
intersections of tropical varieties. One can force two tropical varieties to intersect 
in the correct dimension by translating one of them locally by a generic vector. 
The intersection multiplicities are then computed using the formula from the fan 
displacement rule. 



An intersection product in matroid fans was introduced in S2 , FR . In particular 



this made it now possible to do intersection theory on moduli spaces. 



Tropical intersection theory has many applications. For example, one can use FS2 
to see that certain intersection products in toric varieties can be computed as trop- 
ical intersection products. It has also been used in |BS| to study the relative re- 
alizability of tropical curves. A prominent example of the usefulness of tropical 
intersection theory is enumerative geometry (see for example GKM , Rl , KM| ) 



One can formulate many combinatorially complex enumerative problems in terms 
of much simpler intersection products on moduli spaces. 

However, in concrete cases these products are still tedious to compute by hand, 
especially in higher dimensions. For the purpose of testing new conjectures or 
studying examples of a new and unfamiliar object, one is often interested in such 
explicit computations. This paper aims to analyse how one can efficiently compute 
divisors, products of tropical cycles and other constructions frequently occurring in 
tropical intersection theory. 

After briefly discussing the basic notions of polyhedral and tropical geometry, we 
study some basic operations in tropical geometry in Section [3] We describe how a 
lattice normal vector and the divisor of a rational function are computed and we 
give an algorithm that can check whether a given tropical cycle is irreducible. 

In Section [4] and [5] we discuss algorithms to compute intersection products in K n 
and to create matroid fans. 
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Section [6] is the main focus of this paper: Here we discuss how the combinatorial 
structure of the moduli spaces A4 n of rational n-marked tropical curves can be 
used to efficiently compute these spaces or certain subcycles thereof. We introduce 
Priifer sequences and show how they are in bijection to combinatorial types of 
rational curves. Since Priifer sequences are relatively easy to enumerate, we can 
make use of this to compute (parts of) Ai n - We also show how the combinatorics 
of a rational curve can be retrieved from its metric vector and we study the local 
structure of A4 n . More precisely, we show that locally around each point A4 n is 
the Cartesian product of some K fc and several Mi with i < n. 

In Section [7] we list some open questions and the main features of a-tint, an 
extension for polymake that implements all of the algorithms discussed in this 
paper. We also include some benchmarking tables that show how some algorithms 
compare to one another or react to a change in certain parameters. 
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2. Preliminaries: Polyhedral and tropical geometry 

In this section we establish the basic terms and definitions of polyhedral and tropical 
geometry needed in this paper. For a more thorough introduction to polytopes and 
polyhedra see for example [Zj. 

2.1. Polyhedra and polyhedral complexes. A polyhedron or polyhedral cell in 
V = M. n is a set of the form 

a = {x ■ Ax > b} 

where A e M. mxn ,b e M. m , i.e. it is an intersection of finitely many halfspaces. We 
call a a cone if b = 0. 

Equivalently, any polyhedron a can be described as 

a = conv{pi, . . . ,p k } + M>o?T + ••• + R> r t + L 

where p±, . . . ,pk, T\, . . . , rj e W 1 and L is a linear subspace of W 1 . The first de- 
scription is often called an ^-description of a and the second is a V-description of 
a. 

It is a well known algorithmic problem in convex geometry to compute one of 
these descriptions from the other. In fact, both directions are computationally 
equivalent and there are several well-known convex-hull-algorithms. Most notable 
are the double-description method |MRTT , the reverse search method AF and 



the beneath-and-beyond algorithm (e.g. [G],[e]). Generally speaking, each of these 
algorithms behaves very well in terms of complexity for a certain class of polyhedra, 
but very badly for some other types (see ABS| for a more in-depth discussion of 



this). In general it is very difficult to say beforehand, which algorithm would be 
optimal for a given polyhedron. It is also an open problem, whether there exists 
a convex hull algorithm with polynomial complexity in both input and output. 



All of the above mentioned algorithms are implemented in polymake GJ . At the 
moment, all algorithms in a-tint use the built-in double-description algorithm that 
was implemented by Fukuda [F]. 
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For any polyhedron a we denote by V a the vector space associated to the affine 
space spanned by cr, i.e. 

V<j := (a - b; a, b e cr) 

We denote by A CT := V a n Z™ its associated lattice. The dimension of cr is the 
dimension of V a . 

A face of cr is any subset r that can be written as o~nH, where H = {x : (x, a) = A} is 
an affine hyperplane such that cr is contained in one of the halfspaces {x : (x, a) > A} 
or {x : (x,a) < A} (i.e. we change one or more of the inequalities defining a to an 
equality). If r is a face of cr, we write this r < a (or r < a if the inclusion is proper). 
By convention we will also say that cr is a face of itself. 

Finally, the relative interior of a polyhedron is the set 

rel int (cr) := cr \ (^J t 

T<(7 

A polyhedral complex is a set E of polyhedra that fulfills the following properties: 

• For each a e E and each face r < cr, r e E 

• For each two cr, cr' e E, the intersection is a face of both. 

If all of the polyhedra in E are cones, we call E a fan. 

We will denote by E^ fc ^ the set of all /c-dimensional polyhedra in E and set the 
dimension of E to be the largest dimension of any polyhedron in E. The set- 
theoretic union of all cells in E is denoted by |E|, the support of E. We call E 
pure- dimensional or pure if all inclusion-maximal cells are of the same dimension. 
We call E rational if all polyhedral cells are defined by inequalities Ax > b with 
rational coefficients A. If not explicitly stated otherwise, all complexes and fans in 
this paper will be pure and rational. 

Note that a polyhedral complex is uniquely defined by giving all its top-dimensional 
cells. This is the way in which complexes are also usually handled in polymake: You 
specify the vertices and rays of the complex and then define the top-dimensional 
cells in terms of these rays. Hence we will often identify a polyhedral complex with 
its set of maximal cells. 

The last definition we need is the normal fan of a polytope, i.e. a bounded polyhe- 
dron: Let cr be a polytope, r any face of cr. The normal cone of r in a is 

N Tt<7 := {w € E n : {w,t) = max{(w,3;) \ x e a} for all t e r} 

i.e. the closure of the set of all linear forms which take their maximum on r. These 
sets are in fact cones and the collection of these cones is called the normal fan N a 
of cr. 



2.2. Tropical geometry. Let A be a pure d-dimensional rational polyhedral com- 
plex in W 1 . Let cr 6 and assume r < a is a face of dimension d—1. The primitive 
normal vector of r with respect to cr is defined as follows: By definition there is 
a linear form g e (Z") v such that its minimal locus on a is r. Then there is a 
unique generator of Aa-/A T = Z, denoted by u CT / T , such that g(u a / T ) > 0. One can 
also define this for a polyhedral complex in a vector space A ® R, for a prescribed 



lattice A (see for example [GKM ). If not stated otherwise, we will however always 



consider the standard lattice Z™. 

A tropical cycle (A, w) is a pure rational d-dimensional complex X together with 
a weight function u : X^ -*■ Z such that for all codimension one faces r e X^ 1 ^ 
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it fulfills the balancing condition: 

Y, wfcK/r = € V/V T 

(7>T 

We call X a tropical variety if furthermore all weights are positive. 

Two tropical cycles are considered equivalent if they have the same support and 
there is a finer polyhedral structure on this support that respects both weight 
functions. Hence we will sometimes distinguish between a tropical cycle X and a 
specific polyhedral structure X. 

We also want to define the local picture of X around a given cell: Let r e X be any 
polyhedral cell. Let II : V -» V/V T be the residue morphism. We define 

Star* (r) := {M> • H(a - r); r < o e X} u {0} 

which is a fan in V/V T (on the lattice A/A T ). If we furthermore equip Star^(r) 
with the weight function wstar(R>o ' n(cr - r)) = Uj(a) for all maximal a, then 
(Star*(r), cjgtar) is a tropical fan cycle. 

3. Basic computations in tropical geometry 

3.1. Computing the primitive normal vector. The primitive normal vector 
u a j T defined in the previous section is an essential part of most formulas and cal- 
culations in tropical geometry. Hence we will need an algorithm to compute it. 
An important tool in this computation is the Hermite normal form of an integer 
matrix: 

Definition 3.1. Let M e Z mxn be a matrix with n > m. We say that M is in 
Hermite normal form (HNF) if it is of the form 

M = (0 mx („_ m ),T) 

where T = (tjj) is an upper triangular matrix with i^, > and for j > i we have 

Remark 3.2. We are actually only interested in the fact that T is an upper tri- 
angular invertible matrix. Furthermore, it is known that for any A e Z mx ™ there 
exists a U £ GL„(Z) such that B = AU is in HNF (see for example [C| 2.4]). 

Proposition 3.3. Let X c R n be a d- dimensional tropical cycle, t e X^ 1 ^. Let 
A e ^(» l -< i + 1 ) >< ™ suc h that V T = kerA and V a = kerA, where A denotes A without its 
first row. Let U e GL„(Z) such that 

AU = (0( n _d+i) x (d-i),T') 

is in HNF. Denote by U*i the i-th column of U . Then: 

(1) [7*i, ... , Uid-i is a Z-basis for ker A=V T 

(2) [7*i, ... , U*d is a 2,-basis for ker A = V a 

In particular U*d = ±u a / T mod V T . 

Proof. It is clear that U*i, . . . , U*d-i form an K-basis for kerA and the fact that 
det U = ±1 ensures that it is a Z-basis. Removing the first row of A corresponds to 
removing the first row of AU so we obtain an additional column of zeros. Hence 
[7*1, . . . , U*d is a basis of ker A and U*d is a generator of A CT /A T . □ 
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Remark 3.4. In [Cj 2.4.3], Cohen suggests an algorithm for computing the HNF of 
a matrix using integer Gaussian elimination. However, he already states that this 
algorithm is useless for practical applications, since the coefficients in intermediate 
steps of the computation explode too quickly. A more practical solution is an LLL- 
based normal form algorithm that reduces the coefficients in between elimination 
steps, a-tint uses an implementation based on the algorithm designed by Havas, 
Majevski and Matthews in |HMM| . 

Note that, knowing the primitive normal vector up to sign, it is easy to determine 
its final form, since we know that the linear form defined by u a j T must be positive 
on a. So we only have to compute the scalar product of U*d with any ray in a that 
is not in r and check if it is positive. 

What remains is to compute the matrix A such that A = V T . There is an obvious 
notion of an irredundant "H-description of r: Assume 

r = f]H i nf]S j 

i=l j = l 

where Hi = {x e M. n ; (x,Zi) = a*} and Sj = {x e ]R n ; (x,Wj) > f3j} for some Zi,Wj e 
7j n ,ai,/3j e K. This is considered an irredundant "H-description if we cannot remove 
any of these without changing the intersection and we cannot change any of the 
inequalities into an equality. Note that most convex hull algorithms return such an 
irredundant description. Now it is basic linear algebra to see the following: 

Lemma 3.5. Let r 6 W 1 be a polyhedron given by an irredundant %-representation 
T = f]H i nf]S j with Hi = {ieR"; (x, z 4 ) = aj. Denote by H° = {x e R n ; {x,z,) = 0}. 
Then 

r /V 

K = n#° = H ; 

Furthermore, if r is a codimension one face of a polyhedral complex and r < a, then 
there is an I e {1, . . . , r} such that 

Remark 3.6. There is an additional trick that can make lattice normal computa- 
tions speed up by a factor of up to several hundred. Assume you want to compute 
normal vectors of a 10-dimensional variety in K 20 . In this case we would have to 
compute the HNF of 11 x 20-matrices. However, for computing u a j T we can project 
a onto V a . Now the matrix of the codimension one face r is only a 1 x 10-matrix. 
The normal form of this matrix can of course be computed much faster. Note that 
we have to take care that the projection induces a lattice isomorphism on a. For 
this, we have to compute a lattice basis of u, which still requires computation of 
an HNF of the matrix associated to a - but only once, instead of once for each 
codimension one face. 



3.2. Divisors of rational functions. The most basic operation in tropical inter- 
section theory is the computation of the divisor of a rational function. Let us first 
discuss how we define a rational function and its divisor. Our definition is the same 
as in 



AR 



Definition 3.7. Let X be a tropical variety. A rational function on A" is a function 
ip : X -> E, that is affine linear with integer slope on each cell of some arbitrary 
polyhedral structure X of X. 
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polymake example: Computing a tropical variety. 

This creates the weighted complex consisting of the four orthants of M. 2 with weight 
1 and checks if it is balanced. The maximal cones are represented in terms of indices 
of the rays (starting the count with 0). 

atint > $w = new WeightedComplex( 

RAYS=>[[1,0] , [-1,0] , [0,1] , [0,-1]] , 

MAXIMAL_CDNES=> [ [0 , 2] , [0 , 3] , [1 , 2] , [1 , 3] ] , 

TR0PICAL_WEIGHTS=> [1,1,1,1]); 
atint > print $w->IS_BALANCED; 
1 



The divisor of tp on X, denoted by tp-X, is defined as follows: Choose a polyhedral 
structure X of X such that tp is affine linear on each cell. Let X' = ;f(dimX-i) be 
the codimension one skeleton. For each r e X' , we define its weight via 

^ v -x(r) = I Y, w ( cr )^(" CT /r) J - ¥>T I X) u( a ) u <ylr I 
\<j>t I \<r>T I 

where tp„ and <p T denote the linear part of the restriction of tp to the respective cell. 
Then 

tp-X := (X^uj^x) 

Remark 3.8. While the computation of the weights on the divisor is relatively easy 
to implement, the main problem is computing the appropriate polyhedral structure. 
The most general form of a rational function tp on some cycle X would be given 
by its domain, a polyhedral complex Y with |X| £ \Y\ together with the values and 
slopes of tp on the vertices and rays of Y. To make sure that tp is affine linear on 
each cell of X, we then have to compute the intersection of the complexes, which 
boils down to computing the pairwise intersection of all maximal cones of X and 
Y. Here lies the main problem of computing divisors: One usually computes the 
intersection of two cones by converting them to an "H-description and converting 
the joint description back to a V-description via some convex hull algorithm. But as 
we discussed earlier, so far no convex hull algorithm is known that has polynomial 
runtime for all polyhedra. Also, [T] shows that computing the intersection of two 
V-polyhedra is NP-complete. 

Hence we already see a crucial factor for computing divisors (besides the obvious 
ones: dimension and ambient dimension): The number of maximal cones of the 
tropical cycle and the domain of the rational function. Table [l] in the appendix 
shows how divisor computation is affected by these parameters. 

Example 3.9. The easiest example of a rational function is a tropical polynomial 

tp(x) = max{(ui, x) + a,-; i = 1, . . . ,r} 

with Vi 6 1i n ,cti € K. To this function, we can associate its Newton polytope 

P v = conv{(w 4 , en); i = 1, . . . , r} c R n+1 

Denote by N v its normal fan and define 2V* := N v n {x : x n+ \ = 1}. Then AT* can 
be considered as a complete polyhedral complex in M. n and it is easy to see that 
tp is affine linear on each cone of this complex. In fact, each cone in the normal 
fan consists of those vectors maximizing a certain subset of the linear functions 
(vi, (x\, . . . ,x n )) + on • x n+ \ at the same time. 

So for any tropical polynomial tp and any tropical variety X we can compute an 
appropriate polyhedral structure on X by intersecting it with AT* An example is 
given in Figure [T] 
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Figure 1. The surface is X = max{l, x, y, z, -x, -y, -z} ■ R 3 with 
weights all equal to 1. The curve is max{3a; + 4, x-y-z, y+ z+3}-X , 
the weights are given by the labels. 



polymake example: Computing a divisor. 

This computes the divisors displayed in figure [T] 



atint 


> 


$f = new MinMaxFunction( 






INPUT _STRING=> "max ( 1 , x , y , z , -x , -y , -z) ") ; 


atint 


> 


$x = divisor (linear _nspace (3) ,$f ) ; 


atint 


> 


$g = new MinMaxFunction( 






INPUT _STRING=> "max (3x+4, x-y-z, y+z+3) ") ; 


atint 


> 


$c = divisor ($x,$g) ; 



3.3. Irreducibility of tropical cycles. A property of classical varieties that one is 
often interested in is irreducibility and a decomposition into irreducible components. 
While one can easily define a concept of irreducible tropical cycles, there is in general 
no unique decomposition (see Figure [2]). We can, however, still ask whether a cycle 
is irreducible and what the possible decompositions are. 

Definition 3.10. We call a (i-dimensional tropical cycle X irreducible if any other 
d-dimensional cycle Y with \Y\ £ \X\ is an integer multiple of X. 




Figure 2. The curve on the left is irreducible. The curve on the 
right is reducible and there are several different ways to decompose 
it. 

To compute whether a cycle is irreducible, we have to introduce a few notations: 

Definition 3.11. Let X be a tropical cycle with a fixed polyhedral structure X. 
Let N be the number of maximal cells o"i,...,crjv of X. We identify an integer 
vector u) e Z N with the weight function <7j Wj. We define 
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• Ax '■= e %> N ■ (X,uj) is balanced} (which is a lattice). 

• V x ■= A x ® K 

Now fix a codimension one cell r in X. Let 5 be the induced polyhedral structure 
of Star^(r). For an integer vector to e Z™, we denote by tug the induced weight 
function on S. Then 

• A T X := {oj e 7L N : (S,u>s) is balanced} 

Remark 3.12. We obviously have Ax = n re ^( dimX - 1 ) ^x an d similarly for Vx- 
Clearly, if X is irreducible, then dim Vx should be 1 and vice versa (assuming that 
gcd(wi, . . . , wjv) = 1, where the are the weights on X). However, so far this 
definition is tied to the explicit choice of the polyhedral structure. We would like 
to get rid of this restriction, which we can do using Lemma 3.15| Hence we will 



also write Vx and Ax- We call Vx the weight space and Ax the weight lattice of 
X. 

Definition 3.13. Let (X, oj) be a d-dimensional tropical cycle and X a polyhedral 
structure on X. We define an equivalence relation on the maximal cells of X 
in the following way: Two maximal cells a, a' are equivalent if and only if there 
exists a sequence of maximal cells a = oo, . . . , a r = a' , <jj £ X'^ such that for all 
i = 0,...,r-l, the intersection er^ n cr,j +1 is a codimension one cell of X , whose only 
adjacent maximal cells are Oi and <Ji + \. 

Lemma 3.14. Let a and a' be equivalent. Then: 

(1) Va = Vo< 

(2) uj(a) = uj(<t') 
Proof. 

We can assume that ana' -.re ^(dimX-i) 

(1) Choose any representatives Vo-Jt-jV^^ of the lattice normal vectors. Then 
w(cr)w (T / r +uj(a')v a , /r e V T 
Let gi , . . . , g r e A v such that 

(gA 



V a = ker 



\9r) 



Since V T £V a , we have for all i: 

= g l {u{a)v <J i T + u(a')Vv'i T ,) 
= u}(a')gi(v a ,/ T r) 

Hence v a ,j r e V a and since V a < = V T x (v a ,/ T ), we have V CT < £ V a . The other 
inclusion follows analogously. 
(2) X is balanced at r if and only if Star^(r) is balanced, which is a one- 
dimensional fan with exactly two rays. Such a fan can only be balanced if 
the weights of the two rays are equal. 

□ 

Lemma 3.15. Let X and X' be two polyhedral structures of a tropical cycle X. 
Then Vx = Vx 1 (and similarly for Ax = Vx n 1t N ). 
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Proof. We can assume without loss of generality that X' is a refinement of X . De- 
note by N and N' the number of maximal cells of X and X' , respectively and fix 
an order on the maximal cells of both structures. First of all, assume two maximal 
cones of X' are contained in the same maximal cone of X. Since subdividing a 
polyhedral cell produces equivalent cells in terms of definition |3.13| they must have 



the same weight by Lemma 3.14 Thus the following map is well-defined: We par- 
tition {1, . . . , N'} into sets Si,..., Sn such that j e Si <=>• Oj £ <y% (where erj and 
<7j are maximal cells of X' and X , respectively). Pick representatives . . . , jn} 
from each partitioning set Si and let p : Vx' -* be the projection on these coor- 
dinates jk- By the previous considerations, the map does not depend on the choice 
of representatives. We claim that lm(p) £ Vx- Let r be a codimension one cell of X 
and t' any codimension one cell of X' contained in t. Then Star^(r) = Star^'(r'), 
so if lj e 7L N makes X' balanced around r', then p(uj) makes X balanced around 
r. Bijectivity of p is obvious, so Vx = Vx'- 

□ 

Theorem 3.16. Let(X,uj) be a d- dimensional tropical cycle. Then X is irreducible 
if and only if g ■= gcd(w(cr), a e X^) = 1 and dim Vx = 1. 

Proof. Let X be irreducible. Clearly g must be 1, since otherwise a rational multiple 
of X would provide a full-dimensional cycle in X not equal to k ■ X for an integer 
k. Now assume dimVx = rank(Ax) > 1. Then we have an element uj' e 1 N that is 
not a multiple of to and such that (X,w') is balanced, which is a contradiction to 
our assumption that X is irreducible. 

Now let g = 1 = dimVx- Assume X is not irreducible. Then we can find a poly- 
hedral structure of X and two weight functions co',ui" on this polyhedral structure 
such that (X,uj'),(X,oj") are both balanced and uj' + k ■ u>" for any integer k. In 
particular dimVx > 2, which is a contradiction. □ 

After having laid out these basics, we want to see how we can actually compute 
this weight space: 

Proposition 3.17. Let r be a codimension one cell of a d-dimensional tropical 
cycle X in W 1 . Let ui, . . . Uk £ Z" be representatives of the normal vectors u a j T for 
all a > t. Also, choose a lattice basis l\,...,ld-i of A r . We define the following 
matrix: 

M T := {u x ... u k l x ... l d _ x )z7. n < k+d ~V 

Then k T x = 7r(ker(A/ T ) nZ fe+rf_1 ) x 7LS N ~ k ^ , where n is the projection onto the first 
k coordinates and N is again the number of maximal cells in X . 

Proof. Fix an order on the maximal cells of X and let 

J := {j e [N] : t is not a face of <7j}. 

Then clearly lS N ~ k > = (e^; j e J) z £ A^ and it is easy to see that A^- must be iso- 
morphic to Z( N ~ k *> x A Starx ( T ). Hence it suffices to show that A starx ( r ) is isomorphic 
to 7 r(ker(Af r )nZ fe+rf - 1 ). 

Let (a 1 ,...,a k ,b 1 ,...,k) e ker(A/ T ) n Z fc+d_1 , Then Y,am = T.(-k)k e Ar, so 
Starx(T) is balanced if we assign weights a^. In particular (ai, . . . ,afc) € Ag tarjf ( T ). 
Since l\,...,ld-i are a lattice basis, any choice of the such that Starjsf(T) is 
balanced fixes the bi uniquely, so 7r is injective on ker(M T ) and surjective onto 

A-Star x (T)- □ 
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Algorithm 1 weightSpace(X) 

l: Input: : A pure-dimensional polyhedral complex X 

2: Output: : Its weight space Vx 

3: V X = R W 

4: for r a codimension one face of X do 
5: Compute M T as above 

6: VZ = 7r(ker(M T )) + (ej ■ r is not a face of aj) 
7: V x = V x nV x 

8: end for 

9: return Vx 



polymake example: Checking irreducibility. 

This creates the six-valent curve from example [2] and computes its weight space (as 
row vectors). 

atint > $w = new WeigrrtedComplex( 

RAYS=>[[1,0] , [1,1] , [0,1] , [-1,0] , [-1,-1] , [0,-1]] , 

MAXIMAL_CDNES=> [ [0] , [1] , [2] , [3] , [4] , [5] ] , 

TR0PICAL_WEIGHTS=> [1,1,1,1,1,1]); 
atint > print is_irreducible ($w) ; 

# FALSE is displayed as an empty result 
atint > print cycle_weight_space ($w) ; 
1-110 
1 1 
10 10 
10 10 



This finally allows us to give an algorithm that computes Vx- 

Remark 3.18. One is often interested in the positive weights one can assign to a 
complex X to make it balanced. This is now very easy using polymake: Simply 
intersect Vx with the positive orthant (K^o)^ and you will obtain the weight cone 
of X. 



4. Intersection products in R™ 
There are two main equivalent definitions for an intersection product in K", the 



fan displacement rule [FS2 Theorem 3.2] and via rational functions AR . At first 
sight, the computationally most feasible one seems to be the latter, since we can 
already compute it with the means available to us so far: 

Let X,Y be tropical cycles in W 1 and ipi = max{xi,yi} : R n x M" R. Denote by 
7r : M" x M n M." the projection onto the first n coordinates. Then we define 

X-Y:= 7T*(Vl ^ n -(XxY)) 

(Here, applying 7r* just means forgetting the last n coordinates) However, comput- 
ing this directly turns out to be rather inefficient. The main reason is that, since we 
compute on the product X x Y, we multiply the number of their maximal cones by 
each other and double the ambient dimension. As we have discussed earlier, both 
are factors to which the computation of divisors reacts very sensitively. 

A different definition of the intersection product is given by Jensen and Yu in |JY| : 
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Theorem 4.1. Let X,Y be tropical cycles in R™ of dimension k and I respectively. 
Let a be a (k + I - n) -dimensional cone in the complex X nY and p any point in 
rel int (a). Then a is a cell in X - Y if and only if the Minkowski sum 

Starx(p) - Stary(p) 

is complete, i.e. its support is R™. 



This definition is very close to the fan displacement rule and it is in fact not difficult 
to see that they are equivalent. So, at first glance it would seem to be an unlikely 
candidate for an efficient intersection algorithm. In particular, it is in general 
algorithmically undecidable, whether a given fan is complete (see for example the 
appendix of [n]). However, one can also show that Starx(p) -Stary(p) can be made 



into a tropical fan (see JY for more details). Since R ra is irreducible, a tropical 
fan is complete if and only if it is n-dimensional. In this case it is a multiple of R". 

The weight of the cone a in the above Theorem is then computed in the following 
manner: 

Lemma 4.2 ( |JY ). Let a be a polyhedral cell in X - Y . Let p e rel int (a). Then 
uj X -y(o-)= Y, w x (p 1 )-w y (p 2 )-((A pi +A P2 ) : A pi+P2 ) 

pi eStarx (p) ,P26Stary (p) 
s.t. perel int(pi-p2) 



This now allows us to write down an algorithm based on these ideas: 



Algorithm 2 MinkowskiIntersection 

l: Input: Two tropical cycles X,Y in R n of codimension k and I respectively, 

such that k + I < n 
2: Output: Their intersection product X ■ Y 

3: Compute the (n - (k + Z))-skeleton Z of X n Y 

4: for cr a maximal cell in Z do 

5: Compute an interior point p e rel int a 

6: Compute the local fans Starx(p), Stary(p) 

7: if for any p\ e Starx(p),P2 s Stary(p) the cell p\ - pi is n-dimensional then 

8: Compute weight lux-y of a as described above 

9: else 

10: Remove a 

11: end if 

12: end for 

13: return (Z,ljx-y) 



polymake example: Computing an intersection product. 

This computes the self- intersection of the standard tropical line in R 2 . 

atint > $1 = tropical_lnk(2,l) ; 

atint > $i = intersect ($1, $1) ; 

atint > print $i->TROPICAL_WEIGHTS ; 
1 
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5. Matroid fans 

Matroid fans are an important object of study in tropical geometry, since they are 
the basic building blocks of what we would consider as "smooth" varieties. There 
are several different but equivalent ways of associating a tropical fan to a matroid, 



see for example AK FR FS1 S3 S . One possibility, which immediately implies a 
method to compute the fan, is given in [FS1| Proposition 2.5]: 

Definition 5.1. Let M be a matroid on n elements. For w e W 1 let M w be the 
matroid whose bases are the bases a of M of maximal w-cost £ i6(T Wi. Then w lies 
in the Bergman fan B(M) if and only if M w has no loops, i.e. the union of its bases 
is the complete ground set. 

Remark 5.2. The convex hull of the incidence vectors of the bases of a matroid is a 
polytope in K™, the so-called matroid polytope Pm- So the vectors w maximizing a 
certain basis are exactly the vectors in the normal cone of the vertex corresponding 
to that basis. Hence the Bergman fan is a subfan of the normal fan of Pm- In addi- 
tion, we know that it has dimension rank(M) (this follows immediately from other 
possible definitions of B(M), see for example [ AK| ) . This gives us an algorithm to 
compute B(M): 



Algorithm 3 bergmanFanFromNormalFan 

l: Input: A matroid M on n elements, given in terms of its bases. 

2: Output: Its Bergman fan in 1" 

3: Compute the normal fan F of the matroid polytope Pm- 

4: S = the rank(M)-skeleton of F. 

5: for £ a maximal cone in S do 

6: Let p be the corresponding face of Pm maximized by £ 

7: Let <r\ , . . . , o~d be the bases corresponding to the vertices of p. 

8: if U <Ji £ [n] then 

9: Remove £ from S 
10: end if 

11: end for 

12: return (S, w = 1) 



While this algorithm is fairly simple to implement, it is highly inefficient for two 
reasons: Computing the skeleton of a fan from its maximal cones can be fairly 
expensive, especially if we want to compute a low-dimensional skeleton. But mainly, 
the problem is that from the potentially many cones of S we often only retain a 
small fraction. Hence we compute a lot of superfluous information. 

5.1. Computing matroid fans via circuits. A different definition of a matroid 
fan can be given in terms of its circuits: 

Definition 5.3. Let M be a matroid on n elements. Then B(M) is the set of all 
elements w e K™ such that for all circuits C of M, the minimum min-fu^; % e C} is 
attained at least twice. 

This definition actually produces the inversion in of the fan in the first definition, 



but that is obviously not a problem. It is used by Felipe Rincon in R2 to compute 
the Bergman fans of linear matroids, i.e. matroids associated to matrices. His algo- 
rithm requires the computation of a fundamental circuit C(e, I) for an independent 
set / and some element e £ I such that I u {e} is dependent: 

C(e, J) = {e} u {ie \ {i}) u {e} is independent} 
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It is an advantage of linear matroids that fundamental circuits can be computed 
very efficiently purely in terms of linear algebra. For general matroids it can still 
be computed using brute force. With this modified computation of fundamental 
circuits the algorithm of Rincon can be used to compute Bergman fans of general 
matroids. It turns out that this is still much faster than the normal fan algorithm 
above. Table [4] in the appendix demonstrates this. 



polymake example: Computing matroid fans. 

This computes the Bergman fan of a matrix matroid and of the uniform matroid 
C/3,4. 



atint 


> 


$m = 


new Matrix<Rational>( [[1,-1,0,0] ,[0,0,1,-1]]); 


atint 


> 


$bm = 


= Bergman_f an_linear ($m) ; 


atint 


> 


$u = 


matroid: :unif orm_matroid(3,4) ; 


atint 


> 


$bm2 


= Bergman_f an_matroid($m) ; 



5.2. Intersection products on matroid fans. Intersection products on matroid 



fans have been studied in S2 , FR . Both approaches however are not suitable for 
computation. While the approach in |S2| is more theoretical (except for surfaces, 
where its approach might lead to a feasible algorithm), the description in |FR| might 
seem applicable at first. The authors define rational functions which, applied to 
B(M) xB(M), cut out the diagonal. Hence they can define an intersection product 
similar to |AR| . 

However, these rational functions are defined on a very fine polyhedral structure 



of B(M) induced by its chains of flats. These are very hard to compute ( KBE 



gives an incremental polynomial time algorithm but SI states that already the 
number of hyperplanes can be exponential) and the subdivision computed by the 
algorithm of Rincon is in general much coarser. Also, recall that this approach to 
computing an intersection product already proved to be inefficient in W 1 . 

It remains to be seen whether there might be a more suitable criterion for compu- 
tation of matroid intersection products, maybe similar to Theorem |4.1| 

6. Moduli spaces of rational curves 
6.1. Basic notions. We only present the basic notations and definitions related 



to tropical moduli spaces. For more detailed information, see for example GKM 



Definition 6.1. An n-marked rational tropical curve is a metric graph with n 
unbounded edges, labelled with numbers {l,...,n}, such that all vertices of the 
graph are at least three- valent. We can associate to each such curve C its metric 
vector (d(C)jj*)j<j e rW, where d(C)i_j is the distance between the unbounded 
edges (called leaves) marked i and j determined by the metric on C. 

Define $„ : R™ R®, a » (a, + a,-)i<j- Then 

M„ ■■= {d(C);C n-marked curve} £ r(2)/$„(R™) 

is the moduli space of n-marked rational tropical curves. 

Remark 6.2. It is shown (e.g. in |GKM| ) that M. n is a pure (n-3)-dimensional fan 
and if we assign weight 1 to each maximal cone, it is balanced (though they do not 
use the standard lattice, as we will see below). Points in the interior of the same 
cone correspond to curves with the same combinatorial type, i.e. forgetting their 
metric, they are equal. In particular, maximal cones correspond to curves where 
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each vertex is exactly three-valent. We call this particular polyhedral structure on 
M. n the combinatorial subdivision. 

The lattice for M n under the embedding defined above is generated by the rays 
of the fan. These correspond to curves with exactly one bounded edge. Hence 
each such curve defines a partition or split I\I C on {l,...,n} and we denote the 
resulting ray by vi (note that vj = Vic). Given any rational n-marked curve, each 
bounded edge Ei of length on induces some split = 1, . . . d on the leaves. In the 
moduli space, this curve is then contained in the cone spanned by the vi i and can 
be written as Y.QiVif 

While this description ofM n is very useful to understand the moduli space in terms 
of combinatorics, it is not very suitable for computational purposes. By dividing 
out Im($„), we have to make some choice of projection, which would force us to do 
a lot of tedious (and unnecessary) calculations. Also, the special choice of a lattice 
would make normal vector computations difficult. However, there is a different 
representation of M. n '- It was proven in |AK| and |FR| , that 

M„ = £(^-i)/ ((I,---, 1)) R 

as a tropical variety, where K n ~\ is the matroid of the complete graph on n - 1 
vertices. In paricular, matroid fans are always defined with respect to the standard 
lattice. Dividing out the lineality space ((1, . . . , 1)) of a matroid fan can be done 
without much difficulty, so we will usually want to represent M n internally in ma- 
troid fan coordinates, while the user should still be able to access the combinatorial 
information hidden within. The remaining parts of this section will be dedicated 
to this purpose. 

While the description of M n as a matroid fan automatically gives us a way to 
compute it, it turns out that this is rather inefficient. Furthermore, as soon as we 
want to compute certain subsets of M. n , e.g. Psi-classes, the computations quickly 
become infeasible due to the sheer size of the moduli spaces. Hence we would like 
a method to compute M n (or parts thereof) in some combinatorial manner. The 
main instrument for this task is presented in the following subsection. 

6.2. Prufer sequences. Cayley's Theorem states that the number of spanning 
trees in the complete graph K n on n vertices is n n ~ 2 . One possible proof uses so- 
called Prufer sequences: A Prufer sequence of length n—2 is a sequence (ai, . . . , a„_2) 
with aj € {1, . . . , n) (Repetitions allowed!). One can now give two very simple algo- 
rithms for converting a spanning tree in K n into such a Prufer sequence (Algorithm 
[4j and vice versa (Algorithm [5]) . 



Algorithm 4 PrueferSequenceFromGraph(T) 

l: Input: A spanning tree T in K n 

2: Output: A sequence P = (ai, . . . , a n -2),di e [n] 

3: /': (): 

4: while No. of nodes of T > 2 do 

5: Find the smallest node i that is a leaf and let v be the adjacent node 
6: Remove i from T 
7: P=(P,v) 

8: end while 



It is easy to see that this induces a bijection (see for example AZ Chapter 30]). 
An example for this is given in Figure [3] 
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Algorithm 5 GraphFromPrueferSequence(P) 

l: Input: A sequence P = (a 1; . . . , a n _ 2 ), a* 6 [n] 

2: Output: A spanning tree T in if„ 

3: T •■= graph on n nodes with no edges 

4: V := [n] 

5: while |V| > 2 do 
6: Let i := min V \ P 
7: Let j be the first element of P 
8: Connect nodes i and j in T 
9: Remove i from V and the first element from P 

10: end while 

11: Connect the two nodes left in V 




(5,6,5,6) 



Figure 3. An example for converting a spanning tree on Kq into 
a Priifer sequence and back. The tree can also be considered as 
a 4-marked rational curve with additional labels at the interior 
vertices. 

As one can see from the picture, tropical rational n-markcd curves with d bounded 
edges can also be considered as graphs onn + d+1 vertices: We convert the un- 
bounded leaves into terminal vertices, labelled 1, . . . , n and arbitrarily attach labels 
n + 1, . . . ,n + d + 1 to the other vertices. This will allow us to establish a bisec- 
tion between combinatorial types of rational curves and a certain kind of Priifer 
sequence: 

Definition 6.3. A moduli Priifer sequence of order n and length d is a sequence 
(oi, . . . , a n+ d-i) for some d > 0, n > 3 with aj e {n + 1, . . . , n + d + 1} such that each 
entry occurs at least twice. 

We call such a sequence ordered if after removing all occurrences of an entry but 
the first, the sequence is sorted ascendingly. 

We denote the set of all sequences of order n and length d by V n ,d and the corre- 
sponding ordered sequences by V^d- 

Example 6.4. The sequences (6, 7, 8, 7, 8, 6) and (6, 7, 6, 7, 8, 8) are ordered moduli 
sequences of order 5 and length 2, but the sequence (6, 8, 8, 7, 7, 6) is not ordered. 

Definition 6.5. For fixed n and d we call two sequences p, q e V n ,d equivalent if 
there exists a permutation a e §({n + 1, . . . , n + d+ 1}) such that g>j = a(pi) for all i. 

Remark 6.6. It is easy to see that for fixed n and d the set P^ d forms a system of 
representatives of V n ,d modulo equivalence, i.e. each sequence of order n and length 
d is equivalent to a unique ordered sequence. 

We will need this equivalence relation to solve the following problem: As stated 
above, we want to associate Priifer sequences to rational tropical curves by assigning 
vertex labels n + 1, . . . , n + d+ 1 to all interior vertices. There is no canonical way to 
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do this, so we can associate different sequences to the same curve. But two different 
choices of labellings will then yield two equivalent sequences. 

Proposition 6.7. The set of combinatorial types of n-marked rational tropical 
curves is in bisection to U2=o d • More precisely, the set of all combinatorial 
types of curves with d bounded edges is in bisection to d . 

Proof. The bijection is constructed as follows: Given an n-marked rational curve C 
with d bounded edges, consider the unbounded leaves as vertices, labelled {1, . . . , n}. 
Assign vertex labels {n + 1, . . . , n + d - 1} to the inner vertices. Then compute 
the Priifer sequence P(C) of this graph using Algorithm [4] and take the unique 
equivalent ordered sequence as image of C. 

First of all, we want to see that P(C) e V n ,d- Since C has n + d+ 1 vertices if consid- 
ered as above, the associated Priifer sequence has indeed length n + d- 1. Further- 
more, the n smallest vertex numbers are assigned to the leaves, so they will never 
occur in the Priifer sequence. Hence P(C) has only entries in {n + 1, . . . , n + d + 1}. 
In addition, it is easy to see that each interior vertex v occurs exactly val(i>) - 1 
times (since we remove val(w) - 1 adjacent edges before the vertex becomes itself a 
leaf), i.e. at least twice. 

Injectivity follows from the fact that if two curves induce the same ordered sequence, 
they can only differ by a relabelling of the interior vertices, so the combinatorial 
types are in fact the same. Surjectivity is also clear, since the graph constructed 
from any P € d is obviously a labelling of a rational n-marked curve. □ 

We now present algorithm [6j which, given a moduli sequence, computes the cor- 
responding combinatorial type in terms of its edge splits (it is more or less just a 
slight modification of Algorithm [5]) : 

Theorem 6.8. In the notation of Algorithm^ the procedure generates the set of 
splits Ii,...,Id of the combinatorial type corresponding to P. More precisely: If 
v{i) is the element chosen from V in iteration i e {n + 1, . . . , n + d- 1}, then I n -i is 
the split on the leaves {1, . . . ,n} induced by the edge {v(i),pi}. 

Proof. Let v(i) be the element chosen in iteration i, corresponding to vertex uii in 
the curve. In particular i £ P. This means that i has already occurred val(w;) - 1 
times in the sequence P as P[0]. Hence the node v(i) is already (val(wi) - 1)- 
valent, i.e. connected to nodes qj',j = 1, . . . , val(iUi) - 1. If qj is a leaf (i.e. < n), then 
qj e A v f{\. Otherwise, qj must have been chosen as v(k) in a previous iteration 
k < i. Hence it must already be val(u>/c)-valent. Inductively we see that each node 
"behind" qj is either a leaf or has already full valence. In particular, no further 
edges will be attached to any of these nodes. 

By induction on i, the edge {v (i), qj} (assuming qj is not a leaf) corresponds to the 
split A q .. In particular, A q . has been added to A v . Hence A v is the split induced 
by the edge {v(i),pi}. □ 

Example 6.9. Let us apply algorithm [6] to the following sequence P e V$ 5 (see 
figure [4] for a picture of the corresponding curve) : 

P = (9, 9, 10, 10, 11, 11, 12, 12, 13, 13, 14, 14). 

The algorithm begins by attaching the leaves {1, . . . , 8} to the appropriate vertices, 
i.e. after the first for-loop we have Ag = {1,2}, Aw = {3,4}, An = {5,6},Ai 2 = {7,8} 
and P= (13,13,14,14),!/= {9,. . .,14}. Now the minimal element oiV\Pisv = 9. 
We set I\ = Aq = {1, 2} to be the split of the first edge. Then we connect the vertex 9 
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Algorithm 6 COMBINATORIALTYPEFROMPRUEFERSEQUENCE(P,n) 

l: Input: A moduli sequence P = (pi, . . . ,Pn) 6 V n ,d 

2: Output: The rational tropical n-marked curve associated to P in terms of the 
splits I±,...,Id induced by its bounded edges. 

3: d= N -n+1 

4: V = {l,...,7i + d+l} 

5: A n +i, . . . ,A n+c i+i = 

6: //First: Connect leaves 
7: for i = l. ..n do 
8: A Pz = A p . U {i} 
9: V=V\{i} 
10: P= (p i+1 ,...,p N ) 

11: end for 

12: //Now create edges 

13: for i = n+ 1 . . .n + d- 1 do 

14: v = mmV \ P 

15: 1%—n ~ A v 

16: if length(P) > then 

17: // We denote by P[0] the first element of the sequence P. 
18: A P[Q] = A P[a] u A v 

19: V = V\{i} 

20: P= (pi+i,...,p N ) 

21: end if 
22: end for 

23: //Create final edge 

24: I d = A minV 



to the first vertex in P, which is 13. Hence j4 13 = A^uAg = {1,2}. We remove 9 from 
V and set P to be (13, 14, 14). Now v = minT^ \ P = 10. We obtain the second split 
h = A w = {3,4}. The we connect vertex 10 to 13, so A 13 = A 13 u A w = {1,2,3,4}. 
We set V = {11, 12, 13, 14} and P = (14, 14). In the next two iterations we obtain 
splits ^3 = An = {5,6},/4 = A12 = {7,8} and we connect both 11 and 12 to 14, 
setting A14 = {5, 6, 7, 8}. Now P = () and V = {13, 14}, so we leave the for-loop and 
set the final split to be I5 = A13 = {1,2, 3, 4}. 



3 4 5 6 




Figure 4. The curve corresponding to the moduli sequence P, 
including labels for the interior vertices. 

6.3. Enumerating maximal cones of A4 n . We now want to apply the results of 
the previous section to compute M n - For this it is of course sufficient to compute 
all maximal cones. More precisely, we will only need to compute all combinatorial 
types corresponding to maximal cones, i.e. rational n-marked tropical curves whose 
vertices are all three-valent. Using Algorithm [6] we can then compute its rays 
Vz t , . . . , Vi n _ 3 . These can easily be converted into matroid coordinates with the 
construction given in |FR| Example 7.2]. 
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Proposition 6.7 directly implies the following: 

Corollary 6.10. The maximal cones of M n are in bijection to all ordered Priifer 
sequences of order n and length n — 3, i.e. sequences (oi, . . . , a2n-4,) with ai in 
{n + 1, . . . , In- 2} such that each entry occurs exactly twice. 

This also gives us an easy way to compute the number of maximal cones of Ai n , 
thus reproving the formula of |SS[ Theorem 3.4]: 

Lemma 6.11. The number of maximal cones in the coarse subdivision of Ai n is 

n-4 

Il(2(n-0-5) 

i=0 



Proof. We prove this by constructing ordered Priifer sequences of order n and length 
n — 3. The sequence has 2n - 4 entries. Since it is ordered, the first entry must 
always be n + 1. This entry must occur once more, so we have 2n - 5 possibilities to 
place it in the sequence. Assume we have placed all entries n+\,...,n + k, each of 
them twice. Then the first free entry must be n+k+1, since the sequence is ordered 
and we have 2{n - k) - 5 possibilities to place the remaining one. This implies the 
formula. □ 



As one can see, the complexity of this number is in (D(n n ~ 3 ), so there is no hope 
for a fast algorithm to compute all of M n for larger n (except using symmetries). 
As we will see later, however, we are sometimes only interested in certain subsets 
or local parts of M n - 



polymake example: Computing M. n . 

This computes tropical .Mo, 8 and displays the number of its maximal cones. 

atint > $m = tropical_m0n(8) ; 

atint > print $m->MAXIMAL_CONES->rows () ; 

10395 



6.4. Computing products of Psi-classes. A factor that often occurs in inter- 



section products on moduli spaces are Psi-classes. In KM , the authors describe 
tropical Psi-classes as (multiples of) certain divisors of rational functions, but also 
in combinatorial terms. For nonnegative integers k\, . . . , k n and I £ [n], they define 
K(I) := Y,iei ki- Then one of their main results is the following Theorem 



Theorem 6.12 (KM Theorem 4.1]). The intersection product tp^ 1 tp^-Mn is 

the subfan of M. n consisting of the closure of the cones of dimension n-3-K([n\) 
corresponding to the abstract tropical curves C such that for each vertex V of C we 
have val(l^) = K(Iy) + 3, where 

I v = {i e [n] : leaf i is adjacent to V } £ [n]. 

The weight of the corresponding cone cr(C) is 



UveV(C)K(I V )l 



In combination with Proposition 6.7 this allows us to compute these products in 
terms of Priifer sequences: 
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Corollary 6.13. The maximal cones in f/^ 1 ip 1 ^™ ' M- n o,re in bijection to the 

ordered moduli sequences P 6 „_3_x([n]) ^at fulfill the following condition: 

Let d = n — 3 - X([n]) and ki = for i = n + 1, . . . ,n + d — 1. For any a e 
{n + 1, . . . , n + d + 1} let ji, . . . ,ji( a ) 6 {1, . . . ,n + d - 1} be the indices such that 
Pj i = a. Then 

1(a) 

l(a)=2+£k jt 
i=i 

Proof. Recall that any entry a corresponding to a vertex v a in the curve C(P) 
occurs exactly val(w a ) - 1 times. By the Theorem above the valence of a vertex is 
dictated by the leaves adjacent to it. Furthermore, the leaves adjacent to a vertex 
v a can be read off of the first n entries of the sequence: Leaf i is adjacent to v a if 
and only if P^ = a. 

So, given a curve in the Psi-class product, vertex v a must have valence 3 + K(I Va ), 
so it occurs 2 + K(I Va ) = 2 + X i: p. =Q ^ times. Conversely, given a sequence fulfilling 
the above condition, we obviously obtain a curve with the required valences. □ 

We now want to give an algorithm that computes all of these Priifer sequences. 
As it turns out, this is easier if we require the ki to be in decreasing order, i.e. 
fci > ki > ■ ■ ■ > k n . In the general case we will then have to apply a permutation to 
the ki before computation and to the result afterwards. The general idea is that we 
recursively compute all possible placements of each vertex that fulfill the conditions 
imposed by the ki (if we place vertex a at leaf i with ki > 0, then it has to occur 
more often). Due to its length, the algorithm has been split into several parts: 
iteratePlacements goes through all possible entries of the Priifer sequence re- 
cursively. It uses placements to compute all possible valid distributions of an 
entry, given a certain configuration of free spaces in the Priifer sequence. 



Algorithm 7 psiProductSequencesOrdered^i, . . . , k n ) 
1: Input: Nonnegative integers k\ > &2 > • • • > k n 

2: Output: All Priifer sequences corresponding to maximal cones in ip 1 ^ 1 ip 1 ^ 

M n 
3: K = £ h 

4: current jvertex = n + 1 

5: current .sequence = (0, . . . , 0) € z 2n ~ 4 ~ K 

6: exponents = (k x , . . . , k n , 0, . . . , 0) 6 l? n -^ K 

7-. iteratePlacements (current -vertex, current_sequence) 



Proof, (of Algorithm[7| First of all we prove that placements computes indeed all 
possible subsets J e [m] such that | J| = 2 + T,jej kj. So let J = {a x , . . . , a^} be such 
a set with a\ < ■ ■ ■ < ajv- It is easy to see that in each iteration of the while- loop we 
have \ J\ = i- 1. Let 5 = (2 + Y, j€ j kj) - \J\. 

One can see by induction on 5 that, starting in any iteration of the while loop, 
the algorithm will eventually reach an iteration where i is one smaller. This proves 
termination of placements. 

But we can only reach the iteration where i = if in the previous iteration we have 
tried all indices {1, . . . , m} as first element of J. In particular, there was a previous 
iteration, where we chose I = a\ as first element of J. Now assume we are in the 
first iteration where J = {<ii, . . . , a s }, 1 < s < N. Assuming 8 > 0, we can again only 
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lTERATEPLACEMENTS(current_vertex, current_sequence) 

l: if current _vertex > 2n - 2 - K then 

2: if current -sequence contains no O's then 

3: append current .sequence to result 

4: end if 

5: else 

6: / = {i: current_sequence[i] = 0} 

7: for P 6 PLACEMENTS ( (exp[ j] , i 6 /)) do 

8: v = current .sequence 

9: Place current .vertex in v at positions indicated by P 
10: iTERATEPLACEMENTS(current_vertex+l,w) 
11: end for 
12: end if 



PLACEMENTS(fci, . . . , k m ) 

l: Input: Positive integers fci > ••• > k m 

2: Output: All subsets J £ [to] such that | J\ = 2 + £ J€j kj. 

3: Let J = 0,i = 1 
4: while i > do 
5: if | J| < 2 + kj then 

6: Let I e [to] \ J be minimal such that I > max J and we haven't used it yet 

as z-th element of J. 
7: if There is no such I then 

8: STEPDOWN 

9: else 

10: Add I to the indices we have used as i-th element of J 

11: i = i+ 1 

12: J=Ju{l} 

13: end if 
14: else 

15: Add J to solution 

16: STEPDOWN 

17: end if 
18: end while 



STEPDOWN 

l: Forget all indices we used as i-th element of J. 
2: J = J \ {max(J)} 
3: i = i - 1 



decrease i if we have tried all valid placements, including a s+ i. So assume 5=0. 
Then {ai , . . . , a s } is a valid placement, i.e. s = 2 + Y,t=i k ai ■ If we subtract this from 
the equation for J, we obtain 

AT 
i=s+l 

In particular, since the fc, are ordered, we must have k as+1 > 1 and hence also k aj > 1 
for all j < s. This implies 

s 

s = 2 + Y J k az >2 + s 

i=l 
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which is obviously a contradiction. 

With this it is now easy to see that psiProductSequencesOrdered computes 
indeed all the required sequences. □ 

Example 6.14. We do indeed need that k\ >■•■ >k n to be able to compute all 
sequences. Assume n = 7 and (ki . . . kj) = (0, 0, 0, 0, 0, 1,1). A valid sequence would 
be (7, 7, 8, 8, 9, 7, 7, 9), but this sequence would never occur in the algorithm: After 
having placed the first two 7's in placements we would already have 5 = 0, so the 
last two 7's are never tried out. 

For completeness we also give the algorithm for the general case: 



Algorithm 8 PSlPRODUCTSEQUENCES(fci, . . . , fc„) 
1: Input: A list of nonnegative integers k = ki, . . . , k n 

2: Output: All Priifer sequences corresponding to maximal cones in ip^ 1 tjj; 

M n 

3: Let a 6 S n such that <j(k) is ordered descendingly 
4: 1= psiProductSequencesOrdered (<r(fc)) 

5: return cr _1 (/) (applied elementwise to the first n entries of each sequence) 



polymake example: Computing psi classes. 

This computes tpf ■ ^>§ ' "06 ■ -Mo,9 (which is a point) and displays its multiplicity. 

atint > $p = psi_product(9, new Vector<Int> (3 ,2, ,0 , , 1 , 0, ,0) ) ; 

atint > print $p->TR0PICAL_WEIGHTS ; 

60 



6.5. Computing rational curves from a given metric. In previous sections 
we computed rational curves as elements of the moduli space given by their cor- 
responding bounded edges, i.e. the vi that span the cone containing the curve. 
Usually, we will be given the curves either in the matroid coordinates of the moduli 
space or as a vector in rW, i.e. a metric on the leaves. It is relatively easy to 
convert the matroid coordinates to a metric, but it is not trivial to convert the 
metric to a combinatorial description of the curve, i.e. a list of the splits induced 
by the bounded edges and their lengths. 

The paper [B] describes an algorithm to obtain a tree from a metric d on a given 
set S that fulfills the four-point- condition, i.e. for all x,y,z,t e [n] we have 

d(x, y) + d(z, t) < max{d(a;, z) + d(y, t),d(x, t) + d{y, z)} 



and BG Theorem 2.1] shows that the metrics induced by semi-labelled trees (es- 
sentially: rational n-marked curves) are exactly those which fulfill this condition. 

Note that we can always assume d(x, y) > for x t y by adding an appropriate 
element from Im($„). More precisely, if we have an element d e R^) that is 
equivalent to the metric of a curve modulo Im($„), there is a k e N such that 
d + k • <& n (H&i) is a positive vector fulfilling the four-point-condition. In fact, if 
m = d+ T,&i®n(ei) is the equivalent metric, then d + Y.i^i + |cti|) < I > n(ei) still fulfills 
the four-point-condition, since adding positive multiples of $„(ei) preserves it. 

Algorithm [9] gives a short sketch of the algorithm described in [B] Theorem 2] . As 
input, we provide a metric d. We then obtain a metric tree with leaves L labelled 
{1, . . . , n} such that the metric induced on L is equal to d. This tree corresponds 
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to a rational n-marked curve: Just replace the bounded edges attached to the leaf 
vertices by unbounded edges. It is very easy to modify the algorithm such that it 
also computes the splits of all edges. 



Algorithm 9 TREEFROMMETRic(d) [b1 Theorem 2] 



l: Input: A metric d on the set [n] fulfilling the four-point-condition. 
2: Output: A metric tree T with leaf vertices L labelled {1, . . . ,n} such that the 
induced metric on L equals d. 

3: Let V = {l,...,n} 
4: while \V\ > 3 do 

5: Find ordered triple of distinct elements (p, q, r) from V, such that 

r) + d(q, r) - d(p, q) is maximal 
6: Let t be a new vertex and define its distance to the other vertices by 



7: If d(t, x) = for any x, identify t and x, otherwise add t to V. 
8: Attach p and q to t. Then remove p and q from V 
9: end while 

10: Compute the tree on the remaining vertices using linear algebra. 



polymake example: Converting curve descriptions. 

This takes a ray from .Mo, 6 (hi its moduli coordinates) and displays it in different 
representations. 

atint > $m = tropical_m0n(6) ; 
atint > $r = $m->RAYS->row(0) ; 
atint > print $r; 
-1 -1 -1 0-1-10-10 

atint > $c = rational_curve_f romjnoduli ($r) ; 
atint > print $c; 

(1,2,3,4) # This means that this ray represents W{i 2,3,4} 

atint > print $c->metric_vector ; 

1 1 1 1 1 1 1 1 * Read as d(l, 2), d(l, 3), . . . , d(5, 6) 



6.6. Local bases of M. n . When computing divisors or intersection products on 
moduli spaces A4 n , a major problem is the sheer size of the fans, in the number 
of cones and in the dimension of the ambient space. The number of cones can 
usually be reduced to an acceptable amount, since one often knows that only a 
handful of cells is actually relevant. However, the ambient dimension of A4 n is 
(2) ~ n = " 2 3 " 6 C( n2 )- Convex hull computations and operations in linear algebra 
thus quickly become expensive. We will show, however, that locally at any point 
+ p € M n , the span of Star>i n (p) has a much lower dimension. Hence we can do 
all our computations locally, where we embed parts of M n in a lower-dimensional 
space. Let us make this precise: 

Definition 6.15. Let r be a <i-dimensional cone of M. n . We define 



where U(t) = Uct> t r °l hit (a). It is easy to see that for any t p e M. n and t the 
minimal cone containing p, the span of Star^vi n (p) is exactly V(r). 




d(t,p) = -(d(p,q) - d(p,r) - d(q,r)) 



d(t, x) = d{x,p) - d{t,p) for x + p 



V(t) :=({a>T;<jeM n }) R =(U(T)) 
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We are now interested in rinding a basis for this space V(t), preferably without 
having to do any computations in linear algebra. The idea for this is the following: 
Let C T be the combinatorial type of an abstract curve represented by an interior 
point of t. We want to find a set of rays Vj, all contained in some a > r, that 
generate V(t). Each such ray corresponds to separating edges and leaves at a 
vertex p of C T along a new bounded edge (whose split is of course I\I C ). We will 
see that for a fixed vertex p with valence greater than 3, all the rays separating 
that vertex span a space that has the same ambient dimension as M va \( p )- I n fact, 
it is easy to see that they must be in bijection to the rays of that moduli space. 

Hence the idea for constructing a basis is the following: In addition to the rays 
of t, we choose a basis for the "M. va x( p )" at each higher- valent vertex p. This 



choice is similar to the one in KM Lemma 2.3]. There the authors show that 
Vk •= {vs,\S\ = 2,k i S} is a generating set of the ambient space of M n for any 
k e [n] and it is easy to see that by removing any element it becomes a basis. 

Now fix a vertex p of C T such that s ■= val(p) > 3. Denote by Ii,...,I s the splits 
on [n] induced by the edges and leaves adjacent to p (in particular, some of the Ij 
might only contain one element). We now define 

W p := {v huIi \i,j * l,i * j} 

(This corresponds to the set V\ described above) and 

B p ■= Wp \ {v l2Ul3 }. 

Clearly all the following results also hold if we choose i,j + k for some k > 1 
or remove a different clement in the definition of B p (in particular, because the 
numbering of the Ii is completely arbitrary). To make the proofs more concise, we 
will however stick to this particular choice. We introduce one final notation: For 
\Ij \ = 1, we set vj j ■■= 0. 

Lemma 6.16 (see also [KM| Lemmas 2.4 and 2.7]). 

(1) Let p be a vertex of the generic curve C T and define I\ , . . . , / s , W p as above. 
Then 

viW p \j>l ) 

(2) Let vj be a ray in some a > t and assume it separates some vertex p of C T . 
Define I\, . . . ,I S ,W P as above. Assume without restriction that I\ c I c . 
Then 

v i = Y, v s - ( TO_ 2 )( Y, y ij + v i) - E vsm.odV T . 



Proof. 



(1) We define a = (at) € R n via o< = 1, if i e ii and ai = (s - 3) otherwise. 
Furthermore we define b = (bi) e M. n via 

0, if i is a leaf attached to p 

1 , if z is not a leaf at p and lies in I\ 
(s - 3), if i is not a leaf at p and does not lie in I\. 
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We now prove the following equation (to be considered as an equation in 
R"), where each ray is represented by its metric vector): 



/ 



E v=(s-3) 



\ 



E 



■W/i -<l>n(b) + <Pn(a)- 



We index by all sets T = {ki,k 2 },fa * k 2 . We have 

if fcl,fc 2 € = l,...,s 



E - 



T 



0, 

s - 2, if ki e , k 2 e Ij , j > 1 
2(s-3), if /ci e Ii,k 2 e j > l;i * j. 



We now study the right hand side in four different cases: 

(a) If ki , k 2 e Ji , then both are not leaves at p. Hence the right hand side 
yields + 0- 2 + 2 = 0. 

(b) If fci, k 2 6 > 1, again both are not leaves at p. The right hand side 
now yields + - 2(s - 3) + 2(s - 3) = 0. 

(c) Assume k\ e Ii,k 2 e Ij,i,j > 1 and i + j. If both are not leaves at 
p, we get 2(s - 3) + - 2(s - 3) + 2(s - 3). If only one is a leaf, we 
get (s - 3) + - (s - 3) + 2(s - 3). Finally, if both are leaves, we get 

+ 0- + 2(s-3). So in any of these cases the right hand side agrees 
with the left hand side. 

(d) Assume fci e I\,k 2 e Ij 7 j > 1. If both are not leaves, we get (s - 3) + 

1 - (s - 3) - 1 + (s - 2). The other cases are similar. 

(2) We know that / must be a union of some of the Ij and we assume without 
restriction that / = Uj>fe Ij for some k > 1. Furthermore we define 

m := \{i : l t c I}\ = s -k + l. 

We now prove the following formula (again in M}- 2 '. A similar formula for 
the representation of a ray vj in Vk and a similar proof can be found in 



KM Lemma 2.71 



v i = E v huij -(m- 2)4> n I E e ' I 

i,j>k ) 



(m - 2) j Yj v Ij +«/ - 4>n (E e i I I 
\j>k \i=l I ) 



(6.1) 



To see that the equation holds, let us first compute w. We index K^) by 
all sets T := /c 2 }, fci * k 2 . Then we have 







(e^ 




\3>k 





1, if fci € /, fc 2 ^ / or vice versa 
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Hence 

0, if {fci, £2} £ h for some i > k 

= • 2, if kx 6 I, f*2 i I or vice versa 

{kiM} [2, tfkieli,k 2 elj,i±j;i,j>k 

1 0, if {fci, k 2 } £ ij for some i > fc 

2, otherwise. 




Finally we get 

( w ){fei,fe 2 } 



I -2, if {ki,k 2 } £ for some i > fc 
I 0, otherwise. 

Thus it remains to prove that 

I 2(m - 2), if {hi, k 2 } £ I; for some i > k 
t 0, otherwise. 



( v '- z \k lM } 



For this, let k\ + k 2 e [n]. If T := {^1,^2} £ ii for some i > k, then 
(w/)r = (viiUi 3 )T = for all i,j > k and (4> n (T,iei e/))r = 2. Thus the 
formula holds. Now if k\ e Ii,k 2 e Ij for i + j and i,j > k, we still have 
(vj)j- = 0. Furthermore, there are (m-2) choices for a ray vj. u j' with j' + j 
and (m - 2) choices for a ray Vi'_ u j. with i' * i. For these rays, the T-th 
entry is 1, for all other rays vj' u j' it is 0. Hence 

i j 

= 2(m-2) = (m-2)U„Q>)j . 

< T \ lei l r 

Finally, if k\ 6 I (say fci e 7^), fc 2 i I, then (vj)j- = 1. There are (m - 1) 
choices for a ray vi i j j with j + i. Since (4> n {Y,ui e/))r = 1) we get 

(m-l) = (m-2)U n (Y i ei)) +1. 

\ lei It 

Hence equation |6 . 1 1 holds . 





□ 



Theorem 6.17. Let ve x , ■ ■ ■ , ve t be the rays of r. Then the set 

B T := |J BpU{u El ,...,D Bt } 



val(p)>3 



is a basis for V(t). In particular, the dimension ofV(r) can be calculated as 
dim^(r) =dimr+ £ (f^^) -val(p) ) . 

val(p)>3 

Proof. By Lemma |6.16| these rays generate V^ T : We can write each in some a > r 
in terms of Wp and the bounded edges at the vertex associated to it. The first 
part of the Lemma then yields that we can replace any occurrence of vi 2U i 3 to get 
a representation in B p and the bounded edges. 

To see that the set is linearly independent, we do an induction on n. For n = 4 the 
statement is trivial. For n > 4, assume r is the vertex of M. n . Then B T actually 
agrees with the set Vk \ {vs} for some S and we are done. So let p be a vertex 
of C T that has only one bounded edge attached and denote by i one of the leaves 
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attached to it. It is easy to see that applying the forgetful map ftj to B T , we get 
the set -Bf ti ( r ). If p is three- valent, then the ray corresponding to the bounded edge 
at p is mapped to and all other elements of B T are mapped bijectively onto the 
elements of _B ft .( r ). Since the latter is independent by induction, so is B T . 

If p is higher- valent, only rays from B p might be mapped to or to the same 
element. Hence, if we have a linear relation on the rays in B T , we can assume by 
induction that only the elements in B p have non-trivial coefficients. But these are 
linearly independent as well: Let q be any other vertex with only one bounded edge 
and j any leaf at q. B p is now preserved under the forgetful map itj and hence 
linearly independent by induction. □ 

At the beginning of this section we introduced the notion that the rays resolving a 
certain vertex of a combinatorial type C T "look like -M val ( p )" . The results above 
allow us to make this notion precise: 

Corollary 6.18. Let M be any polyhedral structure of M n (and hence a refinement 
of the combinatorial subdivision). Let t e M be a d-dimensional cell. Let C' T be the 
combinatorial type of a curve represented by a point in the relative interior of r. 
Denote by p%, . . . ,Pk its vertices and by I the number of bounded edges of the curve. 
Then 

Star_M n (r) = R l ~ d x M val{pi) x ••• x M val(pk) 

Proof. First assume M. is the combinatorial subdivision of Ai n . There is an obvious 
map 

: Star Mn (r) M val{pi) x ••• x M val{Pk) , 
defined in the following way: For each vertex pi of C T fix a numbering of the 
adjacent edges and leaves, Ji,...,/^. Now for each vi in some a > r, there is 
a unique i e {1, ...,fc} such that vi separates Pi. Let S £ such that 

I = Ujes Ij- Again, this choice is unique. Now map vj to vs in M va i( Pi ). It is easy 
to see that this map must be bijective. 



First let us see that the map is linear. By Theorem 6.17 we only have to check 



that the map respects the relations given in Lemma |6. 16 But this is clear, since 



analogous equations hold in M n (again, see KM Lemmas 2.4 and 2.7] for details). 

For any set of rays vj 1 , . . . , vj k associated to the same vertex of C T it is easy to see 
that they span a cone in M n if and only if their images do. Now if a > t is any cone, 
we can partition its rays into subsets Sj,j = l,...,m that are associated to the same 
vertex pj . Each of these sets of rays span a cone o~j which is mapped to a cone in 
•Mvai(pj)- Since a = o~\ x • • • x tr m , it is mapped to a cone in A4 va \( Pl ) x • • • x A^ va i( Pm ). 
Hence ip T is an isomorphism. 

Finally, if M. is any polyhedral structure, let r' be the minimal cone of the combi- 
natorial subdivision containing r. Then I = dinn - ' and we have 

Star M „(r) = R l ~ d x Star^fV) 

□ 
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polymake example: Local computations in A4 n . 

This computes a local version of .Mo. 13 around a codimension 2 curve C with a 
single five-valent vertex, i.e. it computes all maximal cones containing the cone 
corresponding to C . a-tint keeps track of the local aspect of this complex, so it 
will actually consider it as balanced. 

atint > $c = new RationalCurve(N_LEAVES=>13, 

INPUT _STRING=>" (2,3) + (2,3,4) + (1,12) + (1,2,3,4,12) + 
(9,10) + (8,9,10) + (11,13) + (8,9,10,11,13)"); 

atint > $m = local_m0n($c) ; 

atint > print $m->MAXIMAL_CONES->rows () ; 

15 

atint > print $m->IS_BALANCED ; 
1 



7. Appendix 



7.1. Open questions and further research. 



7.1.1. More, efficient 'polyhedral computations. As we already discussed in previous 
sections, most polyhedral operations occurring in tropical computations can in gen- 
eral be arbitrarily bad in terms of performance. However, it remains to be seen 
how different convex hull algorithms compare in the case of tropical varieties. So 
far, a-tint only makes use of the double-description method IF]. Also, measure- 
ments indicate that many computations are much faster when all cones involved 
are simplicial (see the discussion in Section 7.3.1). Since we are not bound to a 
fixed polyhedral structure in tropical geometry, it would be interesting to see if 
we gain anything by subdividing all complexes involved before we do any actual 
computations. 



7.1.2. A computable intersection product on matroid fans. A very interesting ques- 
tion is whether one can find another description of intersection products in matroid 
fans that is suitable for computation. One should be able to compute this from 
the bases of the matroid alone (since that is what is usually given). It would also 
be interesting to see if there is some local, purely geometric criterion similar to 
Theorem |4~T1 



7.2. The polymake extension a-tint. All of what has been discussed in the 
previous sections has been implemented by the author in an extension for polymake 
(www.polymake.org). It can be obtained under 

https : / /bitbucket . org/hampe/atint 

Installation instructions and a user manual can be found under the "Wiki" tab. 
We include here a list of most of the features of the software: 

• Creating weighted polyhedral fans/complexes. 

• Basic operations on weighted polyhedral complexes: Cartesian product, k- 
skeleton, amne transformations, computing lattice normals, checking bal- 
ancing condition. 

• Visualization: Display varieties in M 2 or K including (optional) weight 
labels and coordinate labels. 

• Compute degree and recession fan (experimental). 
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• Rational functions: Arbitrary rational functions (given as complexes with 
function values) and tropical polynomials (min and max). 

• Basic linear arithmetic on rational functions: Compute linear combinations 
of functions. 

• Divisor computation: Compute (k-fold) divisor of a rational function on a 
tropical variety. 

• Intersection products: Compute cycle intersections in M. n . 

• Intersection products on matroid fans via rational functions defined by 
chains of flats (this is very slow). 

• Local computations: Compute divisors/intersections locally around a given 
face/ a given point. 

• Functions to create tropical linear spaces. 

• Functions to create matroid fans (using a modified version of TropLi by 
Felipe Rincon |R3| ) . 

• Creation functions for the moduli spaces of rational n-marked curves: Glob- 
ally and locally around a given combinatorial type. 

• Computing with rational curves: Convert metric vectors / moduli space 
elements back and forth to rational curves, do linear arithmetic on rational 
curves. 

• Morphisms: Arbitrary morphisms (given as complexes with values) and 
linear maps 

• Pull-backs: Compute the pull-back of any rational function along any mor- 
phism. 

• Evaluation maps. Compute the evaluation map ev^ on the labelled version 
of the moduli spaces M n . 

7.3. Benchmarks. All measurements were taken on a standard office PC with 8 
GB RAM and 8x2.8 GHz (though no parallelization took place). Time is given in 
seconds. 

7.3.1. Divisor computation. Here we see how the performance of the computation 
of the divisor of a tropical polynomial on a cycle changes if we change different 
parameters. In Table [T] we take a random tropical polynomial / with I terms, 
where I e {5, 10, 15}. We compute the divisor of this polynomial on X, where X is 
L x M fe_1 in W 1 , L the standard tropical line, i.e. X has 3 cones. We do this ten 
times and take the average. 

Note that the normal fan of the Newton polytope of / usually has around 15 
maximal cones, but is highly non-simplicial: It can have several hundred rays. If, 
instead of /, we take a polynomial whose Newton polytope is the hypercube in R™ 
(here the normal fan is simplicial), then all these computations take less then a 
second. This is probably due to the fact that convex hull algorithms are very fast 
on "nice" polyhedra. 

7.3.2. Intersection products. We want to see how computation of an intersection 
product compares to divisor computation. If we apply several rational functions 

/i, ...,//. to R n , we can compute /i fk • R n in two ways: Either as successive 

divisors /i • ( f 2 •(•■•• K™ ) ) or as an intersection product ( /i • W 1 ) (f k -R n ). Since 

successive divisors of rational functions appear in many formulas and constructions, 
it is interesting to see which method is faster. Table [2] compares this for k = 2. 
We take / and g to be random tropical polynomials with 5 terms and average 
over 50 runs. As we can see, the intersection product is significantly faster in low 
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Table 1 . Divisor of a random tropical polynomial 



dimensions, but its computation time grows much more quickly: For n = 8, the 
intersection product takes seven times as long as the divisors. 



n 


/•(<?•»") 


(/•• ; ." )•(;/•••-•.") 


3 


0.62 


0.14 


4 


0.68 


0.24 


5 


1.04 


0.38 


6 


1.42 


0.84 


7 


1.5 


2.7 


8 


1.6 


11.66 



Table 2. Comparing successive divisors to intersection products 



7.3.3. Matroid fan computation. Here we compare the computation of matroid fans 
with different algorithms. In Table [3] we compare computation of the moduli space 
Ai n , first as the Bergman fan B(K n -i) of the complete graph on n— 1 vertices 
using the TropLi algorithms, then combinatorially as described in Corollary |6.10| 

Table [4] shows some more examples of Bergman fans. We compare the performance 
of the TropLi algorithms R2 to the normal fan Algorithm[3] First, we compute the 
Bergman fan of the uniform matroid U Uj k- Note that we compute it as a Bergman 
fan of a matroid without making use of the matrix structure behind it (the uniform 
matroid is actually realizable). Then we compute two linear matroids, i.e. we let 
TropLi make use of linear algebra to compute fundamental circuits. Ci has as 
column vectors the vertices of the i-dimensional unit cube. 



We did not use the original TropLi program but an implementation of the algorithms in 
polymake-C++. In the case of linear matroids the original program is actually much faster. This 
is probably due to the fact that the data types in polymake are larger and that the linear algebra 
library used by TropLi is more efficient. 
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Table 3. Computation of M 
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Table 4. Computation of several matroid fans 



References 

M. Aigner and G.M. Ziegler, Proofs from THE BOOK, Springer, 1998. 

L. Allcrmann and J. Rau, First steps in tropical intersection theory, Math. Z. 264 (2010), 

no. 3, 633-670, available at |arxiv:0709.3705v3| 

F. Ardila and C.J. Klivans, The Bergman complex of a matroid and phylogenetic trees, 
J. Comb. Theory, Ser. B 96 (2006), 38-49, available at |arxiv:math/0311370v2| 

D. Avis, D. Bremner, and R. Seidel, How good are convex hull algorithms, Computational 
Geometry: Theory and Applications 7 (1997), 265-302. 

D. Avis and K. Fukuda, A Pivoting Algorithm for Convex Hulls and Vertex Enumeration 
of Arrangements and Polyhedra, Discrete and Computational Geometry 8 (1992), 295- 
313. 

J. -P. Barthelemy and A. Guenoche, Trees and Proximity Representations, Wiley- 
Intcrscicncc, Chichester, 1991. 

E. Brugalle and K. Shaw, Obstructions to approximating tropical curves in surfaces via 
intersection theory (2001), available at arxiv:1110.0533v2 

P. Buneman, A note on the metric properties of trees, Journal of combinatorial theory 
17 (1974), 48-50. 

H. Cohen, A course in computational algebraic number theory, 4th ed., Springer Verlag, 
Berlin, 2000. 

H. Edclsbrunner, Algorithms in combinatorial geometry, Springer- Verlag, 1987. 

G. Frangois and J. Rau, The diagonal of tropical matroid varieties and cycle intersec- 
tions, available at arxiv:1012.3260vl 

K. Fukuda, cdd,cddplus and eddlib homepage, 2002. available at http://www.ifor.math. 
ethz . ch/-f ukuda/cdd_home/ 

E. Feichtner and B. Sturmfels, Matroid polytopes, nested sets and Bergman fans, Port. 
Math. (N.S.) 62 (2005), 437-468, available at |arxiv:math/0411260"l 

W. Fulton and B. Sturmfels, Intersection theory on toric varieties, Topology 36 (1997), 
no. 2, available at arxiv:9403002 

A. Gathmann, M. Kerber, and H. Markwig, Tropical fans and the moduli spaces of 
tropical curves, Compos. Math. 145 (2009), no. 1, 173-195, available at |arxiv:070 8.2268 
E. Gawrilow and M. Joswig, polymake: a Framework for Analyzing Convex Polytopes, 
Polytopes — Combinatorics and Computation, 2000, pp. 43-74. polymake is available at 
http : //www . polymake . org 

B. Griinbaum, Convex polytopes, 2nd ed., Springer- Verlag, 2003. 

G. Havas, B.S. Majewski, and K.R. Matthews, Extended gcd and Hermite normal form 
algorithms via lattice basis reduction, Experimental Mathematics 7 (1998), 125-136. 
A.N. Jensen and J. Yu, Stable intersections of tropical varieties, work in progress. 
M. Kerber and H. Markwig, Intersecting Psi-classes on tropical A4o,m Int. Math. Res. 
Notices 2009 (2009), no. 2, 221-240, available at |arxiv:070 9.3953v2 



A- TINT: ALGORITHMIC TROPICAL INTERSECTION THEORY 



31 



[KBE + ] L. Khachiyan, E. Boros, K. Elbassioni, V. Gurvich, and K. Makino, On the complexity 
of some enumeration problems for matroids, SIAM J. Discrete Math. 19 (2006), no. 4, 
966-984. 

[N] A. Nabutovsky, Einstein structures: Existence versus uniqueness, Geometric and Func- 
tional Analysis 5 (1995), no. x, 76-91. 
[M] Grigory Mikhalkin, Tropical geometry and its applications, Proceedings of the ICM, 

Madrid, Spain (2006), 827-852, available at |arxiv:math/0601041v2| 
[MRTT] T.S. Motzkin, H. Raiffa, G.L. Thompson, and R.M. Thrall, The double description 

method, Contributions to theory of games, Vol. 2, 1953. 
[Rl] J. Rau, Intersections on tropical moduli spaces (2008), available at arxiv:0812.3678vl 
[R2] F. Rincon, Computing Tropical Linear Spaces, Journal of Symbolic Computation, to 

appear, available at arxiv:1109.4130 
[R3] , TropLi: a software tool for computing tropical linear spaces. |http: //math. | 

berkeley.edu/~felipe/tropli/ 
[SI] P.D. Seymour, A note on hyperplane generation, Journal of Combinatorial Theory, Series 

B 1 (1994), no. 1, 88-91. 
[S2] K. Shaw, A tropical intersection product in matroidal fans (2010), available at 

[arxiv:1010.3967vT1 

[S3] D. Speyer, Tropical linear spaces, SIAM J. Discrete Math. 22 (2008), 1527-1558, avail- 
able at arxivima th/0410455 
[SS] D. Speyer and B. Sturmfels, The Tropical Grassmannian, Advances in Geometry 4 
(2004), no. 3, 389 - 411, available at math/0304218. 
[S] B. Sturmfels, Solving systems of polynomial equations, CBMS Regional Conferences 
Series in Mathematics, vol. 97, Published for the Conference Board of the Mathematical 
Sciences, Washington, DC (2002). 
[T] H.R. Tiwary, On the Hardness of Computing Intersection, Union and Minkowski Sum 

of Polytopes, Discrete & Computational Geometry 40 (2008), no. 3, 469-479. 
[Z] G. Ziegler, Lectures on polytopes, Springer- Verlag, 1994. Graduate Texts in Mathematics 
152. 

Simon Hampe, Fachbereich Mathematik, Technische Universitat Kaiserslautern, Post- 
fach 3049, 67653 Kaiserslautern, Germany 



E-mail address: hampeOmathematik.uni-kl.de 



